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There have been growing evidences for the critical roles of metal-coordination complexes in defining 
structural and mechanical properties of unmineralized biological materials, including hardness, toughness, 
and abrasion resistance. Their dynamic (e.g. pH -responsive, self-healable, reversible) properties inspire 
promising applications of synthetic materials following this concept. However, mechanics of these 
coordination crosslinks, which lays the ground for predictive and rational material design, has not yet been 
well addressed. Here we present a first-principles study of representative coordination complexes between 
metals and catechols. The results show that these crosslinks offer stiffness and strength near a covalent bond, 
which strongly depend on the coordination state and type of metals. This dependence is discussed by 
analyzing the nature of bonding between metals and catechols. The responsive mechanics of 
metal-coordination is further mapped from the single-molecule level to a networked material. The results 
presented here provide fundamental understanding and principles for material selection in 
metal-coordination-based applications. 



Many biological materials with outstanding mechanical, adhesive and frictional properties benefit from 
their hierarchical composite structures, within which the reinforcing role of metal complexation in the 
polymer network phase was recently highlighted 1 ' 2 . For example, it was reported that the Fe-catecholate 
(1,2-dihydroxybezene) complexes in mussel byssus cuticles provide effective crosslinking strategies between 
metal ions and catecholic amino acid dopa (3,4-dihydroxypheylalanine) in mussel foot protein 1 (mfp-1) 
chains 1 ' 3 . Granules with high density of these crosslinks resemble the particle-reinforced phase in the composite, 
and has been correlated to the hardness and high extensibility of the mussel byssus threads. The metal-dopa 
crosslinks in the matrix were also proposed to endow self-healing properties. Measurements of the mechanical 
properties after removing metals from the cuticle showed 2-fold decrease in the hardness, suggesting a significant 
impact of the coordination complex on the material properties 3 . In the fangs of the wandering spider Cupiennius 
salei, Zn-His coordination complexes were also found in the matrix of chitin and Zn-containing proteins 2 . 
Connection was then identified between the mechanical properties and Zn/Cl enrichment in this fiber- reinforced 
biocomposite. Moreover, the inhomogeneous metal distribution in the fang might have additional implications in 
its graded mechanical properties, with an over two-fold increase in hardness and 1.5-fold increase in stiffness 
measured from the base to the tip 2 . 

The aforementioned examples give prominence to the key roles of metal coordination in defining mechanical 
properties of crosslinked biopolymers. Coordination complexes form crosslinks when two or more ligands each 
donate a nonbonding electron pair to empty orbitals in a transition metal ion. A hexadentate mononuclear triple 
(tris-)Fe-dopa coordination complex was proposed for the mfp-1 crosslinks in the byssus coating 3 , while bi-, 
quadr- and hexa-dentate complex can be induced in synthetic materials by tuning pH 4 . In comparison with other 
crosslinking strategies in polymeric networks, such as non-bonding (e.g. van der Waals, electrostatic, hydrogen 
bonds) 5 and covalent bonds 6 , metal-coordination-based crosslinks play a unique role as they are strong (com- 
parable to covalent bonds), responsive to the environment (e.g. pH), and yet reversible (to prevent catastrophic 
material failure) 4 ' 7 " 9 . These features have been proposed to offer biological materials a number of desirable 
mechanical properties in the absence of mineralization, including enhanced hardness, adhesion, toughness, 
tunability, self-healing capability, efficient energy dissipation upon impact loading, and functions as sacrificial 
bonds 10 " 12 . 
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Figure 1 | Biological and bio-inspired materials with mechanical properties defined by crosslinked biopolymers. (a) Mechanical properties of 
biological materials (e.g. mussel byssus 1 ' 4 , spider fang 2 ) and bioinspired synthetic polymer networks (Fe- and Ti-gels 4 ) are attributed to the presence of 
unmineralized metals that form crosslinking complexes between biopolymers with funcntional groups such as catechol and dopa. (b) Illustration of a 
polymeric network crosslinked by metal-coordination complex consisting of metals and catechols, (c) Atomic structures of (from left to right) covalent 
bonds, bis-catecholate- and tris-catecholate-complexes with coordinate metal-oxygen bonds between metal and catechols. (Images of the mussel byssus 
curticle and spider fang were taken by Genevieve Anderson and Ken Walker, respectively.) 



Recently efforts have been not only made to understand biological 
materials, but also to develop synthetic materials following these 
design strategies (Fig. 1). Possible modifications of the proteins 
and their interfaces were pursued to establish demanding mechanical 
properties 13 . Metal- coordination hydrogels were successfully synthe- 
sized by several groups 41415 , which display distinctly dynamic (e.g. 
responsive, self-healable, tunable, etc.) material properties in con- 
trast to conventional polymers. One excellent example is the pH- 
responsive gel inspired by the pH jump that proteins experience 
during maturation of a mussel byssus secretion, which features high 
stiffness approaching covalently crosslinked gels with additional self- 
healing capability 414 . In order to develop and optimize bio-inspired 
materials following the recipes found in nature and explore the limit 
of maximizing material performance, the understanding of metal- 
coordination-based crosslinks, especially their mechanical signifi- 
cance, is urged 16 . A few key questions can be asked readily: (1) what 
is the nature of the high strength featured by metal-coordination 
complexes that form in specific chemical environment (e.g. pH)? 
(2) to what extend the selection of metals (and ligands) would modu- 
late the mechanical properties of the complex? and (3) how would the 
mechanics of single complex be mapped to the material (e.g. a cross- 
linked polymer network) level? 

To address these issues, we present here first a molecule-level 
study of the coordination complexes. First-principles calculations 
were performed to investigate representative coordinate and cova- 
lent crosslinks, with focus on their structural and mechanical prop- 
erties. The impact of these crosslinks to mechanics of materials at the 
network level is then explored, which reveals the responsive hallmark 
of metal-coordination based biological and bio-inspired materials. 



Results 

We considered both metal-catecholate coordination complexes and 
covalent bonds as possible crosslinking mechanisms in biopolymer 
networks. Ultraviolet-visible (UV) absorption spectra studies 
showed that at low pH (<7), catechol is reduced and coordination 
complexes are dissolved instead of crosslinked 17 . At higher pH, dou- 
ble (bis-) and triple (tris-) complexes form 17 . On the other hand, 
covalent bonds could form by an oxidative process under appropri- 
ate (pH 8.5) conditions 14 . Our first-principles calculations start from 
atomic structures illustrated in Fig. lc. In the bis-complexes, a planar 
square coordination complex was assumed, while the metal ion 
locates in the center of an octahedron for tris-complexes 18 . Five 
metals, including Fe, Mn, Ti, Zn, and Ca, were included in our 
calculations according to their wide appearance in biological (Fe, 
Mn, Zn, Ca) 1 ' 219 and synthetic (Fe, Ti) 41415 polymer materials. Our 
geometry optimization results show that all metals form distorted 
planar square and octahedron configurations with their oxygen 
donors provided by surrounding catechols 18 . The lengths of metal- 
oxygen (MO) bonds / M o were found to be 1.85 (Fe-O), 1.88 (Mn-O), 
1.92 (Ti-O), 2.03 (Zn-O), 2.01 (Ca-O) A in the bis-complexes, and 
1.88 (Fe-O), 1.90 (Mn-O), 2.00 (Ti-O), 2.13 (Zn-O), 2.38 (Ca-O) A in 
the tris-complexes, respectively. The O-M-O bond angles are dis- 
torted from 90° (the angle in undistorted planar squares and octa- 
hedrons) for both bis- and tris-complexes (see Table TI and T2 in the 
Supporting Information for listed values of bond lengths and 
angles). In general, these values are within the range of data reported 
for their oxides 20 . For example, the Ti-O bond lengths obtained here 
(1.92 and 2.0 A) are comparable to the value for Ti0 2 (1.9-2.0 A) and 
Ti 2 0 3 (2.0-2.1 A) 21 . The covalent bond thus forming between two 
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catechol groups was also explored, suggesting a length of 1.49 A for 
the inter- catechol C-C bond. 

The coordinate bonds are commonly considered to be comparable 
to, but weaker than covalent bonds. When the polymer network are 
loaded mechanically as mimicked in this work by pulling apart cate- 
chol groups, the metal ion-catechol crosslinks bear tensile loads. The 
coordination complex thus not only critically defines the stiffness, 
strength, viscoelasticity, and other overall mechanical behaviors, but 
also being responsible for the reversibility and tunability of the net- 
worked material. Here we quantified the mechanics of crosslinks by 
displacing and constraining selected carbon atoms in catechols (see 
Methods section for more details). The force carried by the complex 



was calculated from the reaction force on the loaded or fixed atoms, 
and the results are plotted as a function of relative displacement 
(defined as change in the distance between pulled and fixed carbon 
atoms) in Fig. 2. It was found that in comparison to the covalent 
crosslink, the mechanical response of coordination complexes fea- 
ture a rather complicate profile with multiple peaks, arising form a 
larger number of degree of freedom involved in the deformation 
(bond stretching, bond angle bending, etc.). For bis-complex, the 
ridge was first flattened with negligible responses in the force ampli- 
tude. Further stretching led to subsequent breaking of MO bonds (see 
Fig. 3 and Movie Ml and M2 in Supplementary Information). The 
irregular profiles of force-displacement curves in Fig. 2, especially 
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Figure 2 | Single-molecule mechanics of crosslinks from first-principles calculations. Panels (a-e) plot force- displacement relationships for metal 
(Fe, Mn, Ti, Zn, Ca) -catechol complexes, and (f) for the covalent crosslink. Metal-coordination shows tensile strength of the complex on the order of a few 
nN, varying by the type of metals. 
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Figure 3 | Structural distortion and failure of coordination complexes. We plot here changes in the metal-oxygen bond lengths (labeled in the insets) in 
the metal-coordination complexes with respect to the loading distance (a, c: Ti, b, d: Fe). 



those for tris-complexes, can be explained from the dynamics of MO 
bond distortion and detailed breaking sequence. The Fe- catechol bis- 
complex first experienced a structural distortion, corresponding to 
the first peak, and then failed by breaking two Fe-O bonds simulta- 
neously, while for Ti-catechol bis-complex, Ti-O bonds on the pulled 




Figure 4 | Depedence of crosslinking strength on coordination state and 
metal types. Tensile strength of coordination complexes based on 
different metals (Fe, Mn, Ti, Zn, Ca) and catechols are plotted and 
compared with a covalent crosslink between two catechol groups. It should 
noted that the stength of complexes are contributed by two metal- oxygen 
bonds and thus the metal coordination could reach up to half of the 
strength in a covalent bond. 



catechols broke first. The tris-complexes also showed different bond 
breaking patterns under tensile loading. For Fe-catechol tris-com- 
plex, the two pairs of metal-catechol crosslinks at the left side of the 
metal ion (Fig. lc) broke in two steps, resulting in failure at a much 
larger distance than the Ti-catechol tris-complex, where the breakage 
of crosslinks right to the metal ion led directly to structural failure 
(Fig. 3). These distinct behaviors suggest different stabilities of 
metal-coordination complexes at nonequilibrium conditions, 
although their structures at equilibrium are quite similar as described 
above. 

To further assess the roles of metal ions and their coordination 
states, the whole spectrum of tensile strength for metal- coordination 
complexes was calculated and summarized in Fig. 4. For bis -com- 
plexes, the strength defined as the peak force carried by the crosslink 
consisting of two MO bonds ranges widely from 3.30 (Ca-O), to 8.18 
(Ti-O) nN. Notably, the Fe-catecholate bis-complex features a tensile 
strength of 6.22 nN, while Ti- coordinate crosslink has a tensile 
strength of 8. 18 nN. It should be noted here that this strength defined 
for the coordination complex is actually contributed by a pair of M-O 
bonds that break almost simultaneously (Fig. 3). As a result, single 
MO bond bears a strength approximately half of that value, i.e. — 3.11 
and —4.09 nN for Fe and Ti respectively. This value is close to half of 
the strength of a covalent C-C bond (7.79 nN) as calculated here, but 
several orders higher than the strength of weak non-bonding inter- 
actions such as hydrogen bonds, which are usually in the pN 
range 22 ' 23 . For the tris-complex under tension, the M0 6 octahedron 
was firstly distorted strongly, being elongated along the loading dir- 
ection, then dissociated after the metal ion- oxygen bonds break, 
leaving a bis-complex structure behind (Movie M3 and M4). This 
dissociation defines the strengths of complexes, ranging from 1.89 
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(Ca-O) to 4.87 (Ti-O) nN. The dependence of tensile strength on 
both the coordination state and metal ion types shares the same trend 
as bis-complexes, i.e. Ti-based coordination complex are the stron- 
gest, followed by Fe and Mn featuring similar performance, while Zn 
and Ca show lower values of the strength. Moreover, we could extract 
the tensile stiffness from the linearly elastic part of the curve (Fig. SI). 
The results show that for bis-complexes, Fe-based and Ti-based 
complex could offer stiffness k of 75.61 and 78.28 N/m respectively, 
also close to the covalent crosslink with k = 75.77 N/m, which is one 
order higher than the hydrogen bonds 24 . Similarly, the stiffness of 
single M-O bond is half of the stiffness measured for the coordination 
complex. 

The structural and mechanical properties of the coordination 
complexes can be clarified further from their electronic structures. 
Bader atomic charge analysis 25 was performed to inspect charge dis- 
tribution in the complexes. This approach divides the grid-based 
charge density from quantum mechanical calculations of electrons 
into atomic volumes. The analysis shows that for bis-complexes, 1.3, 
1.7, 1.9, 1.3, and 1.4 electrons are transferred from the metals (Fe, 
Mn, Ti, Zn, Ca) to the catechol groups. For tris-complexes, the 
amount of electrons transferred is 1.3, 1.4, 1.8, 1.3, and 1.5 electrons 
respectively. These results suggest more ionic nature for divalent 
metal such as Ca, which is correlated to the mechanical strength as 
identified above. These results elucidate the combined electrostatic 
and covalent nature of the metal-coordination. To quantify the dif- 
ference in the metal -catechol bonding characteristics, topological 
analysis based on electron localization functions (ELFs) was per- 
formed for all complexes. ELF is a localized function of the 
ground- state electron density and wavefunction obtained from the 
first-principles calculations 26 . The value of ELF ranges from 0 to 1, 
where 1 corresponds to the perfect localization as in covalent bonds, 
and 0.5 corresponds to the electron-gas-like pair probability as in 
metallic bonds. The bonding nature between metal ions and oxygen 
atoms in catechols depends on the type of metal. As shown in Fig. 5a 
and 5b, the results illustrate the nature of coordination bonding 
compared to the covalent and ionic bonds 26 . Moreover, there is a 
distinct contrast in the bonding nature between Ti-O and Ca-O 
bonds. Further comparative analysis shows that only Ca in the five 
metals explored here displays distinctly the ionic nature in the coor- 
dination complex. 

The superior mechanical performance revealed above for single 
metal-catecholate complex could inspire synthetic materials by 
incorporating metal coordination. However, how could these merits 
of individual crosslinks be mapped to the material level remains as a 
critical question to answer. Especially, the contrastive stiffness of the 




Figure 5 | Constractive nature of of metal (Ti, Ca)-catecholate bonding in 
the coordination complex. Isosurfaces at ELF = 0.5 are plotted here with 
the atomic structures of tris-complexes. The distribution of ELF shows 
distinctly covalent and ionic nature for (a) titanium and (b) calcium 
respectively, which explains their difference in the mechanical properties 
(Fig. 2). 



complex to non-bonding interactions opens promising opportun- 
ities in designing highly responsive materials. We now discuss this 
point using a crosslinked polymer network mode with the assump- 
tion of affine deformation. The mechanical response of a typical 
crosslinked polymeric network can be described using a few repres- 
entative units (Fig. 6a). When the representative units are loaded, 
both the polymer chain and crosslink bears an internal force/. We 
define the effective tensile stiffness of polymers (considering their 
intrinsic elasticity, thermal fluctuation, and network entanglement) 
to be K, and the stiffness of crosslinks to be k respectively. At its 
mechanical equilibration, the nature of coordinate bonds could be 
modified into a new state with stiffness k\ e.g by tuning the pH, type 
or concentration of metal ions. Then a response in the force ampli- 
tude is generated inside the material. Straightforward analysis shows 
that response is/ =fa(l + b)l{a + b), where a = k' Ik is the stiffen- 
ing factor of crosslinks, and b = K/k is the relative stiffness between 
the polymer and crosslink. Typical scaling behavior of /// with 
dependence on the stiffness enhancement b is plotted in Fig. 6b. 
The results show that significant response can be established in a 
stressed network, which increases with both a and b. For b = 10, a 
ten-fold stiffening of the crosslink (a = 10, e.g. from hydrogen bond- 
ing to metal- coordination) could yield nearly ten-fold response in the 
internal force carried by the material. We also find that the response 
is more amplified in the softening situation, i.e. with a < 1. 

Discussion 

It was shown in the results that the highest strength achieved in Ti- 
catecholate complex could reach half of the covalent strength. 
Recently, atomic force microscopy (AFM) measurements suggest 
that the strength of a covalent bond is —2.0 nN (for silicon-carbon) 
and —1.4 nN (gold-sulfur) 23 . The mechanics of metal coordination 
in biological materials was characterized as well in AFM 7 and surface 
force apparatus (SFA) 9 studies. The results show that a single dopa 
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Figure 6 | Responsive mechanics of crosslinked networks, (a) Typical 
crosslinking structures in a stressed polymer network can be captured in a 
simple spring model, loaded by force / (b) Response in the internal forces 
(f If) carried by a pre-stressed networks (the inset), with dependence on the 
change in tensile stiffness {k' Ik) of coordinate crosslinks. Kis the effective 
stiffness of the polymer, and k is the stiffness of the crosslink. 
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residue contacting a wet Ti0 2 surface reveals high strength 
(—0.8 nN) yet fully reversible, noncovalent interaction 7 . Moreover 
Fe ions can mediate the unusually strong interaction between posi- 
tively charged proteins through metal coordination with dopa 9 . 
Although there are some theoretical and experimental evidences 
showing that mono- and bis-dentate coordination could form 
between catechol and a Ti0 2 surface, as a result of dissociative 
adsorption 27 , the detailed single-molecule-level chemistry at inter- 
faces in the AFM and SFA experiments is not known. Moreover, the 
interfacial interaction and stability of metal-coordination complex 
could be perturbed by its chemistry and environment such as the 
presence of water 28 . In spite of these uncertainties that prevent rig- 
orously quantitative comparison, our results that are on the same 
order as the reported data remain to be conclusive. 

To further explain the quantitative discrepancy between our cal- 
culations and experimental measurements, it should be noted that 
the failure of metal coordination complexes has dissociation rates 
over a range of time scales. Both thermal fluctuation and loading rate 
will strongly affect the bond dissociation. These effects can be essen- 
tially captured by the Bell-Evans model f h = k B T\n[rx h /(k 0 k B T)]/x h , 
where k B is the Boltzmann factor, T is the temperature, r is the 
loading rate, x h is the equilibrium bond length, and k 0 is the equilib- 
rium dissociation rate 7 . Recent experiments have shown the evidence 
that the stiffness and strength of synthetic gels containing metal- 
coordination crosslinks are rather a function of the dissociation rate 
than the equilibrium constant 15 ' 29 . Our quasi-static loading process 
thus approaches the limit of r = 0 and T = 0. In the force microscopy 
tests, although finite loading rate would yield a higher strength than 
that obtained in our calculations, the thermal fluctuation is expected 
to play a more critical role, which significantly reduces the lifetime of 
metal-ligand bonds and thus the measured strength of coordination 
complexes. This may explain the relatively low values reported in 
experiments. Additionally, although specific metal-coordination 
complex was modeled here according to experimental evidence 17 , 
explicit effects from the chemical environment (pH, ionic strength, 
salt, etc.) could be critical. For example, water or hydroxyl could fill in 
the empty octahedral sites in the bis-complex, ions and molecules 
around could induce charge transfer and polarization as well. The 
solvent effect for metal ions in aqueous solution could vary for dif- 
ferent metals, and thus lead to further contrastive mechanical prop- 
erties of the complexes 30 . However, due to the high computational 
demands to include these effects in first-principles calculations, we 
leave them to be exploited in our future work. 

The nature of coordinate bonds is predominantly electrostatic and 
partially covalent, which differs for different combination of metal 
and donor atoms. There are some understandings from a survey of 
metals in protein structures 31 , that the covalent contributions follow 
a certain pattern - K + < Na + < Ca 2+ < Mn 2+ < Fe 2+ < Co 2+ ~ Ni 2+ 
~ Cu 2+ > Zn 2+ . In comparison to the metal -donor bond lengths 
between metal and oxygen in main- chain carbonyl reported in Ref. 31 , 
our results for Zn and Ca are in good agreement (2.13 vs 2.07 A for 
Zn and 2.38 vs 2.36 A for Ca in tris- complex), while values of Fe-O 
and Mn-O bond lengths are shorter (1.88 vs 2.04 A for Fe and 2.09 vs 
2.19 A for Ca in tris -complex). This may arise because the transition 
metal ions with d n configurations are very sensitive to the perturba- 
tion in their chemical environment that was not include in this work. 
Although rigorous comparison between our first-principles calcula- 
tions results and these referred experimental data cannot be made, 
the covalent bond length range for all complexes and the more dis- 
tinct ionic nature identified for Ca- coordination are in reasonable 
agreement. 

In addition to the catechol considered in this work as the electron 
donor, other orf/zo-dihydroxy moieties (e.g. dopa, titron) 17 and spe- 
cific amino acids (e.g. His, Ser, Thr, Asn, Gin) could also coordinate 
with metals 17 ' 31 . For example, recent work showed that Histidine-rich 
proteins serve as crosslinking points with nitrogen in the spider 



fang 2 . Moreover, experimental evidences were reported for poly- 
chaete jaw 19 ' 32 , Alzheimer Ap aggregation 33 , marine and micro- 
organisms 34 ' 35 , where other transition metal-coordination strategies 
play important roles as well. By comparing and understanding these 
crosslinking mechanisms, extended studies focusing on the material 
selection (in terms of metal and ligands) are expected to be important 
in many practical aspects including material science and biomedical 
applications 36 " 38 . 

Finally at the material level, copolymers containing catechol units 
were synthesized recently, afforded a dual cross-linking system based 
on completely different coordination and oxidation chemistries, and 
competing crosslinking mechanisms 14 . It was shown that oxidation 
yields a stable covalently crosslinked gel after excessive acetic acid 
was added, while metal coordination leads to reversible transition 
between crosslinking and decrosslinking states. Our affme network 
model suggests the tunability of network mechanics by controlling 
either the crosslink (k, k') or polymer (K). In experiments, change 
from k to k' could be established by modulate the metal- coordination 
state, and the effective stiffness K could be elevated from entropy- 
controlled value ~ k B T/L 2 for a flexible polymer chain (L is its size), to 
much higher values in semiflexible nanostructures, e.g carbon nano- 
tubes and graphene 16 ' 39 " 41 . By integrating these thoughts, high-per- 
formance responsive materials could be accomplished. 

In summary, we studied the mechanical properties of metal-cate- 
cholate coordination that play critical roles in defining the perform- 
ance of unmineralized biological materials. Near-covalent stiffness 
and strength of metal-coordination complexes were found, which 
depends on the coordination state and type of metals. The mech- 
anical response in polymers crosslinked by metal coordination was 
predicted through an affine network model. The results, in combina- 
tion with reported experimental evidences that the chemical envir- 
onment can tune the coordination state, implying pronounced 
dynamic mechanical properties. Synthetic materials following this 
concept could thus offer robustness because of their outstanding 
mechanical properties, responsiveness because the metal-coordina- 
tion crosslinks are pH-dependent, reversibility because they can 
break and reform, and tunability because their strength can be modu- 
lated by the metal and chemical environment in question. 

Methods 

Atomic structures of metal-catecholate-coordination complexes. Metal-based bis- 
catecholate- and tris-catecholate-coordination complexes were constructed based on 
the atomic structures of catechols and typical coordination complexes in metal oxides 
(planar square for the bis-complex and octahedron for the tris-complex). Metal- 
catecholate complexes were chosen as our model systems here, because catechol and 
its derivatives are widely present in biological materials, and contained in the side 
chains of dopa 17 . Also catechol serves as a minimal model to reduce the computation. 
Supercells with dimensions of 20, 16, and 16 A were built. The thickness of vacuum 
layers was verified to be large enough to model the complexes as isolated molecular 
systems, by comparing with results using larger supercells. 

First-principles calculations. The structure and properties of the coordination 
complexes were explored using plane-wave-basis-set-based density functional theory 
(DFT) methods. Spin-polarization was included in the calculations. Generalized 
gradient approximation (GGA) in Perdew-Burke-Ernzerhof (PBE) 42 parameter 
settings was used for the exchange-correlation functional and Projector Augmented 
Wave (PAW) 43 potentials were used for ion-electron interaction. We used the Vienna 
Ab-inito Simulation Package (VASP) 44 for the calculations. For all results presented, 
an energy cut-off of 500 eV was used for the plane-wave basis sets. Gamma point was 
used for the Brillouin zone integration. These settings were verified to achieve a total 
energy convergence less than 1 meV per atom. For geometry relaxation, the force on 
atoms converged below 0.01 eVA" 1 . All structures were optimized using the 
conjugated gradient method. This first-principles approach was successfully applied 
to similar systems, e.g. adhesion of catechol on Ti0 2 27 and silica 45 . In previous studies, 
adhesion energy close to experimental measurement was predicted 7 . However, it 
should also be noted that the high computational demand prevents DFT calculations 
to include explicit solvent, and thus the effects of chemical environment (e.g. pH, 
ionic strength, salt, etc.) could not be directly addressed. In this study, this 
information was fed into the molecular models instead, as suggested by absorption 
spectra studies of metal-catecholate complexes 17 . Optimizing the models in first- 
principles calculations could eventually solve this issue by minimizing computational 
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demands while keeping the essential chemistry. This is the subject of our ongoing 
investigation. 

Mechanical tensile tests. Tensile loads were applied on selected carbon atoms in the 
catechols. For the bis-catecholate- and tris-catecholate-complexes, the rightmost two 
carbon atoms (referring to the metal ion, see Fig. lc) were displaced further from the 
other catechols step by step (0. 1 A for each step), while the leftmost two carbon atoms 
(one in each catechol group) were fixed in the loading direction while not constrained 
in other directions. Coordinates of all atoms in the system were scaled accordingly 
after each loading step, to minimize the dependence on the step size. Geometrical 
optimization calculations were performed after each movement, and the force was 
evaluated from the reaction force on displaced or constrained atoms. For the covalent 
bonds between two catechols, the same loading procedure was followed except that 
displacement and constraints were only applied to one atom in each catechol group. 
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